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a Abstract 

g 
I The Stroh formalism is applied to the analysis of infinitesimal surface 

'^ wave propagation in a statically, finitely and homogeneously deformed iso- 

Q tropic half-space. The free surface is assumed to coincide with one of the 

O principal planes of the primary strain, but a propagating surface wave is 

not restricted to a principal direction. A variant of Taziev's technique [Sov. 
Phys. Acoust. 35 (1989) 535] is used to obtain an explicit expression of the 
secular equation for the surface wave speed, which possesses no restrictions 

(T^ on the form of the strain energy function. Albeit powerful, this method does 

^^ not produce a unique solution and additional checks are necessary. How- 

ever, a class of materials is presented for which an exact secular equation 
for the surface wave speed can be formulated. This class includes the well- 

^^ known Mooney-Rivlin model. The main results are illustrated with several 

numerical examples. 
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1 Introduction 

The study of sinall-amplitude surface waves propagating in finitely and homoge- 
neously deformed hyperelastic materials has quite a long history, dating back to 
the classical paper by Hayes and Rivlin [T] in 1961. The interest in using the the- 
ory of small motions superimposed on a large static deformation of a hyperelastic 
half-space is manifold, for once the problem is solved the results are applicable to 
various advanced topics. These, in particular, include the non-destructive evalu- 
ation of solids (see Guz' and Makhort |2| for a review), the incremental stability 
analysis of the loaded surface of a deformed material (see Guz' [3] for comprehen- 
sive review and bibliography), and the acousto-elastic effect (second-order theory 
of linear elasticity) where the pre-deformation is also considered small (see Pao et 
al. ^ for a review of both experimental and theoretical results). 

The results of Hayes and Rivlin, valid for a compressible hyperelastic material, 
have since been extended to the case of hyperelastic materials subject to incom- 
pressibility [S] or to a generic isotropic internal constraint [B] . However, these works 
are limited to the consideration of principal surface waves, i.e. surface waves that 
propagate and attenuate along principal directions of pre-strain. We note that 
Connor and Ogden |7j and Destrade and Ogden |8| considered two-partial (non- 
principal) surface waves polarized in a principal plane of pre-strain. Nevertheless, 
for three-partial non-principal surface waves very few explicit results exist and the 
scope of their applicability is limited. Specifically, Flavin [9] considered the problem 
for a Mooney-Rivlin material with one material parameter much smaller than the 
other; Willson [lOl [11] studied materials subject to equi-biaxial pre-deformations; 
Gerard [12], Gerhart [13], and Iwashimizu and Kobori [13] worked with the lin- 
earized theory of second-order acousto-elasticity; Chadwick and Jarvis [TJ] , Mase 
and Johnson [16] , and Chadwick [17] used the Stroh-Barnett-Lothe integral formal- 
ism; Rogerson and Sandiford [18] used numerical methods for an implicit secular 
equation; etc. 

This paper presents an explicit secular equation for the speed of a surface wave 
propagating in a principal plane, but not in a principal direction, of a tri-axially 
deformed, general hyperelastic incompressible material. This result is achieved by 
using methods first developed by Taziev [191 EO] for surface waves propagating in 
the symmetry plane of a crystal. Although similar, the analysis presented here 
reveals some features particular to the context of nonlinear elasticity. In general, 
the surface wave consists of a linear combination of three partial modes, each 
proportional to exp {ik{n ■ x + qoTri ■ x — vt)}, a = 1, 2, 3, where k is the wave 
number, v the wave speed, n the direction of propagation, m the normal to the 
surface, and g^ an attenuation coefficient. For a general strain energy function, 
the Qa are the roots of a bi-cubic equation with positive imaginary parts. Their 



explicit analytical expressions are awkward and the method of Taziev proves useful 
(Section 4), because it does not require such expressions. However, for a whole class 
of hyperelastic incompressible materials, inclusive of the Mooney-Rivlin model, the 
coefficients of attenuation g^ are obtained analytically, for one is always equal to 
i = a/— T, and so the two others are roots of a bi-quadratic (Section 3). This 
property has been touched upon by Flavin but has been only recently examined 
in detail (Pichugin [Zl])- Here, it leads to the derivation of an exact and explicit 
secular equation for surface waves, which possesses no more than one root for the 
speed. In the general case, the method of Taziev leads to a rationalized secular 
equation (a polynomial of degree 12 in the squared wave speed), with spurious 
roots to be discarded. 

Before these two main results are developed, we summarize in Section 2 the 
basic governing equations and present the equations of motion as a first-order linear 
differential system. The derivation of this system is a lengthy process and is not 
obvious at all. Thanks to some shorthand notations and to the Stroh formalism, the 
system can however be presented in quite a compact form. Its resolution, coupled to 
the appropriate boundary conditions for surface waves (vanishing of the wave away 
from the interface; no incremental traction on the interface), leads to an implicit 
secular equation, to be made explicit in the subsequent two Sections. Section 3 is 
devoted to the special class of incompressible materials which is associated with 
a factorized propagation condition. Numerical examples in this section involve a 
Mooney-Rivlin material characterized by the same material parameters, state of tri- 
axial pre-strain, and normal load as the one considered by Rogerson and Sandiford 
[T8j|. A connection is made with their numerical results for the surface wave speed 
versus the angle of propagation. Some new features are highlighted; in addition, 
the attenuation coefficients for the partial displacements are presented. Section 4 
covers the derivation of the secular equation for the general strain energy density of 
an incompressible material, using a variant of Taziev's technique. To illustrate the 
method, we investigate numerically the case of a deformed Varga material where a 
surface wave propagates in any direction in the plane of shear, with a view to the 
non-destructive acoustic evaluation of a deformed rubber insulator. 

2 Preliminaries 

2.1 Equations of motion 

Consider a half-space, composed of a homogeneous pre-stressed hyperelastic in- 
compressible material with mass density p, characterized by strain energy function 
W. Let (O, Xi, X2, X3) be a fixed rectangular Cartesian coordinate system such that 



the body occupies the region X2 ^ and that the principal stretches coincide with 
the Oxi directions, with corresponding stretch ratios Ai, A2, A3 (Ai 7^ A2 7^ A3 7^ Ai 
and A1A2A3 = 1). The half-space is maintained in the static state of deformation 
by the application of the constant tractions o"i, (T2, a^ at infinity, given by 

a, = XiWi-p, (2.1) 

(no summation assumed here) where Wi := dW/dXi and p is a constant scalar, 
introduced by the constraint of incompressibility. 

Then consider the superposition of a small-amplitude motion u{xi,X2iX^,t) 
upon the primary large deformation. The corresponding incremental nominal stress 
s has the following components 



Sij = BijkiUi^k + puij - p*6ij, (2.2) 

where p* is a Lagrange multiplier, corresponding to an increment in p, and the 
non-zero components of the fourth-order elasticity tensor B are (no summation on 
repeated i,j indexes assumed here) 



iijj ^i^j^* ij^ 



Bi,ij = {Xm - XjW,)Xt/{Xt - A^), 

The linearized incremental equations of motion and incompressibility condition 
read 

Sij,i = puj^tt, Ujj = 0. (2.4) 

We specialize our analysis to the consideration of a surface (Rayleigh) wave, 
propagating in a principal plane but not in a principal direction, see Figure 1. 
Specifically, it is assumed that the wave is traveling with phase velocity v and 
wave number k over the surface X2 = in a direction which makes an angle 6 with 
Oxi, it decays exponentially away from the boundary X2 = 0, and produces no 
incremental traction at the boundary. 

Hence we model this motion by 

{u,p*,s} = {Uikx2),ikP{kx2),ikS{kx2)}e"'^'o'''^'<''''-'"\ (2.5) 

where Ce := cos^, Sg := sin 6', and U, P, S are functions of kx2 alone. By 
substituting these forms of the mechanical displacements and the tractions into the 



incremental constitutive equation (2.2) and using 



of motion and the incompressibility constraint (2.4) can be cast as a homogeneous 



2.1), the incremental equations 



linear system of six first-order differential equations, 

^' = iN^, where ^(A;X2) := [f/i, f/2, t/3, ^21, ^22, ^23]^, (2.6) 




Figure 1: A surface wave propagating in a principal plane but not in a principal 
direction. 

within which the prime denotes differentiation with respect to the variable kx2- 

Here the 6x6 real matrix N follows the usual block decomposition [23] of 
linear anisotropic elasticity, 



N 



N^ N2 

N3 + XI N^ 



X := pv^ 



(2.7) 



where Ni, N2 = N2 5 ^^^ -^3 = -^3 ^^^ 3x3 matrices. However, the components 
of the Ni are specific to the theory of small motions superposed on large static 
deformations [TTj. To present them in a compact form, we introduce the short-hand 
notations 



7,, := {X,W, - X,W,)X^/{X^ - A^) = 7,, + X,W, - X,W, 



J' 



2A,- := X'^iWii - 2XiXjWi, + X]Wjj + 2{X.^VJ - XjW.;)XiXj/{X^, - X]) = 2/3ji. (2.8) 



Then —Ni and N2 are given by 

00(721 - o-2)/72i 

Co se 

59(723 - Cr2)/723 



1/721 


1/723 



respectively, and 



N, 



7] 





— K, 





u 





-K 





fi _ 



(2.9) 



(2.10) 



where 



2ce(/3i2 + 721 - 0-2) + S0731, 

Ce[7l2 - (721 - 0"2)V72l] + Se[732 - (723 
Cg7i3 + 2sl{(323 + 723 - (^2), 

cesei(3i3 - /3i2 - /323 - 721 - 723 + 2^2), 



(^2)^/1: 



23j 



(2.11) 



(see Destrade [21] for the case o"2 = 0, and Destrade and Ogden [S] for the case of 
a wave polarized in the symmetry plane of a stretched and sheared material). 

2.2 Propagation condition 

Now the requirement of exponential decay away from the surface X2 = is ex- 
pressed by choosing the following form for the wave, 

^{kx2) = ee'"''', %q)>0, (2.12) 

where $,° is a constant vector and q is an attenuation coefficient. Then the equa- 



tions of motion (2.6) become the eigenvalue problem: N^" = q^°. The associated 



characteristic equation det(N — ql) = 0, is the propagation condition. This equa- 
tion is a cubic in g^, as demonstrated by Rogerson and Sandiford [T8] : 

72i723g'-[(72i + 723)X-Ci]g4 + (X2-C2X + C3)g2 + (X-C4)(X-C5)=0, (2.13) 

with 



Cl 
C2 
C3 

C4 
C5 



(721713 + 2/3l2723)c0 + (723731 + 2/32372l)s0, 

(723 + 713 + 2/3i2)ce + (721 + 731 + 2/323)s^, 

(712723 + 2/3i27l3)c0 + (721732 + 2/32373i)Se 

+ [712721 + 713731 + 723732 - (/3l3 - /3l2 - /323)^ + 4:(3i2(323\cjsl, 

7i2Ce + 732Se, 



713^ + 2/3l3C0Sg + 73i4- 

Note that the roots qf, g|, gg of the bicubic are such that 



(2.14) 



2 I 2 I 2 

?1 + ^2 + ?3 

2 2 I 2 2 I 2 2 
?ig2 + 92?3 + ?3?1 

2 2 2 
?l92% = 



[(721 + 723)^ - Cl]/(72l723), 

(X2-C2X + C3)/ (721723), 
(X-C4)(X-C5)/(721723). 



(2.15) 



2.3 Implicit secular equation for surface waves 



For the three roots gi, g2, Qs of the propagation condition (2.13) with positive 



imaginary parts, the corresponding eigenvalue problems yield three linearly inde- 
pendent eigenvectors ^^, ^^, ^^, respectively. Then the general solution to the 



equations of motion (2.6) may be written as 



^ikX2) 



(2.16) 



for some constants 71, 72, 73. Explicitly, ^* are given by columns of the matrix 
adjoint to A?" — g^l. Taking, for example, the second such column gives ^* in the 
form: 



a^qf + a2qf + ao 

-Qf + hqf + biQi 
d^qf + d2qf + do 



hqf + hiqi 
{u-X){qf + mqf + n) 

gsQf + giQi 

where expressions for the constants a^, bi, di, hi, and gi are too lengthy to be 
reproduced here and 



(2.17) 



m 



1 (721 - ^2)^ 

721 liii'^-x) 



cl][v-X] + 



+ 



(723 -(^2? 



-2k 



(7 



21 



723 -fisil^-X)' 
0-2)(723-Cr2) 



[/x-X] 



n 



721723 («^ - X) 
(721 - 0^2)^ 2 , (723 - 0^2)^ 2 



-C0S0, 



721 723 

x[(/x-X)(77-X)-«:2]/(72i723) 



(2.18) 



The boundary condition of zero incremental tractions at the plane surface X2 = 
means that 

^(0) = [f/i(0),t/2(0),[/3(0),0,0,0]T. (2.19) 



By comparing this expression with (2.16) evaluated for X2 = 0, we conclude that 



7i) 72; 73 are solutions to a homogeneous linear system of three equations. The 
corresponding secular equation is given by 



[u-X) 



hqf + hqi hsq^ + hiq2 hq^ + hiq^ 
qf + mqf + n gg + "^^i + ^ ^3 + "^?3 + "^ 

93^1 + giQi gsq'i + 9i(i2 g^ql + ^1^3 



0. 



(2.20) 



As noted by Taziev [T^j in the context of linear anisotropic elasticity, this deter- 
minant factorizes greatly. By omitting the factors u — X, qi — qj and hig^ — h^gi, 
we are left with 

nui = ijJiii{m — uii) , (2-21) 

where 

'^i := -(^1 + '?2 + g3), (^u ■= qiq2 + q2q3 + qsqi, ^m ■= -qiq2qi,- (2.22) 



Eq. (2.21) is the secular equation for surface waves in deformed incompressible 
materials; it remains implicit as long as the explicit expressions for the uja are not 
known. 



3 Factorization of the propagation condition 

Before moving on to the derivation of an explicit secular equation, a special case 
must be treated separately. A simple analysis shows that, for a certain class of 
incompressible hyperelastic materials maintained in a state of large homogeneous 



deformation (strain- induced anisotropy), the propagation condition (2.13) factor- 
izes into the product of a term linear in q^ and a term quadratic in g^, thus leading 
to simple explicit expressions for the q^ and eventually, for the secular equation. 
This class of materials includes the well-known Mooney-Rivlin model, often used 
to describe finite deformations of rubber. The described factorization does not in 
general occur for linear elastic materials such as crystals (intrinsic anisotropy). 

3.1 Conditions on the strain energy function and explicit 
secular equation for surface waves 



Following Pichugin pTj, we seek a solution to the propagation condition (2.13) 
of the form q^ = C, where C is a constant independent of v. By substituting 



q = C into (2.15), we obtain three equations for the two quantities S and P, the 



respective sum and product of the remaining q^, namely, 



S + C = [(721 + 723)^ - ci]/ (721723), 

CS + P=iX^- C2X + C3)/(72l723), (3.1) 

CP=-(X-C4)(X-C5)/(72l723). 



After solving (3.1 )i for S and then (3.1)3 for P, the substitutions of S and P into 
(3.1)2 allow for the identification of like-powers of X and Cg on both sides of the 



resulting equation. Since the identity must hold for all X and 6 we obtain that 
C = —1, provided the following relationships are satisfied 



7ij + 7ji = 2l3ij. 
Then the associated propagation condition factorizes into 



(3.2) 



(3.3) 



where 



1 1 \ -- /712 , 7l3\ 2 fl31 , 732 

, + Cg- + 

721 723/ V721 723/ \721 723 

.2 ,, „2\/v ,, „2 ,, 2 



S= ( — + — IX 

P ={X - -fi2ci - 732S0)(X - 7i3Ce - 731 Sg)/ (721723) 



(3.4) 



To sum up, if the conditions ( 3.2[ ) on the strain energy function are satisfied, 
then the bi-cubic (2.13) factorizes into the product of a term linear in q^ and a 



term quadratic in q . It follows that the roots q^ with positive imaginary parts and 



hence the corresponding wi, wn, and um defined in (2.22), can be found explicitly: 
?i = i, '?2g3 = -V?, '?2 + '?3 = iv 2yP - S, so that 



UJl 



i(l + V2VP-^), 



-Uii 



\^+\j2\^-S, uui = iVP. (3.5) 



The end result is that equation ( |2.21 ) is now an explicit secular equation for surface 
waves in deformed incompressible materials satisfying (3.2), namely 



n(l + \/2y/P -S) + ^fp(m + ^P + ^ 2\fP - S 



0. 



(3.6) 



Note that the conditions (3.2) impose restrictions upon the strain energy func- 



tion of the incompressible material. Explicitly, they read {i 7^ j) 

{X} + ?>\])\W,-{?>\} + \])\,W,-{\]-\]){\^,Wu-2\\W,, + X)W,,) = 0. (3.7) 

For example, the Mooney-Rivlin strain energy function, 

W = Vi{Xl + A^ + A^ - 3)/2 + V^iXlXl + A^A^ + XJXJ - 3)/2, (3.8) 

where Vi and V2 are constants, satisfies this condition. This case is treated in the 
following subsection. 

Note also that Pichugin |21j finds conditions under which the propagation con- 



dition (2.13) admits roots of the form q = CX + D, where C, D are constants. It 



turns out this possibility arises when the half-space is subject to an equi-hiaxial pre- 
deformation, whatever the strain energy function may be. Another simplification 



occurs when the bi-quadratic in (3.3) admits a double root (then S = 4P); such is 



the case for the neo-Hookean form of strain energy function (P2 = in (|3.8)). The 



secular equations for surface waves in bi-axially deformed generic incompressible 
materials and in tri-axially deformed neo-Hookean materials were established by 
Willson [1^ and Flavin [9], respectively, and are not investigated further here. 

3.2 Example: Mooney-Rivlin materials 

Now the case of Mooney-Rivlin materials is dealt with in such a way that a con- 
nection is made with the numerical results of Rogerson and Sandiford [18]. For the 



Mooney-Rivlin strain energy function (3.8) the quantities 7jj, /3jj, defined in (2.8), 
yield the following forms, 

7i,- = {Vr + I?2ADA^ 2A, = iV, + V,\l){X^, + A|) = 7ii + 7,., (3.9) 

where k 7^ i, j. 

Using these expressions, together with the material parameters [18] Vi = 2, 
V2 = 0.8, the stretch ratios given by Xj = 3.695, A^ = 0.7, A^ = 0.387, and the 



normal load (J2 = 0.8, the secular equation (3.6) is solved numerically to give the 
variation of the surface wave speed with 9, and the graph of Rogerson and Sandiford 
[18] is reproduced with little effort. We take this opportunity to comment on their 
statement that "the surface wave degenerates into a shear wave as 6 approaches 
and 7r/2" . 

Firstly, we find that at, and close to, the direction Oxi (also the direction of 
greatest stretch), the surface wave speed v = \/X/p is distinct from the bulk shear 
wave \/c4/p, i.e. C5 > C4 in the neighborhood of 6* = 0. For instance, at 6^ = the 
principal surface wave propagation speed a/Xq/p, where y/Xo = 2.917, is found 
from Dowaikh and Ogden's [5] formula Xq = 712 — 72iC^) where 

C^ + C^ + '^^" + ^" ' "'^ ~^'\ - ^^" " "^^' = , (3.10) 

721 721 



(or equivalently, is found from (3.6)), while the bulk shear waves propagate at 
speeds ycs/p and \/ c^j p where ^Jc^ = 1/713 = 3.076 and ^/c^ = ^712 = 2.921. 
Figure 2 shows the variations of these speeds in the (0° — 20°) range. The top 
(dashed) curve is the graph of y/c^, the middle (dotted) curve is the graph of y/cl, 
and the bottom (solid) curve is the graph of y/X. 

Secondly, we find that, as the direction of propagation approaches the Ox^ di- 
rection (direction of least stretch, 9 = 90° here), the surface wave speed v = y/X/p 

10 



^v 




Figure 2: Scaled speeds of bulk and surface waves near the direction of greatest 
stretch in a Mooney-Rivlin material. 



is indeed tending to the bulk shear wave speed \/c^/ p, so that the correspond- 
ing graphs are indistinguishable one from another in the approximative range 
(82° — 90°). However, at 6' = 90° exactly, there exists a two-partial principal surface 
wave whose speed is intermediate between the bulk shear wave speeds \fc^J~p and 
\/c/i/ p. Numerically, when 9 = 90°, ^/c^ = 0.995, a/q = 1.384, and this two-partial 
principal surface wave propagates with speed \/X<^qI p say, where a/Xqo = 1.327, 
the value found from Dowaikh and Ogden's [S] formula Xgg = 732 — 723C^5 iii which 



>3 ^ ^2 ^ 2(/323 + 723 



C^2) - 732 ^ (723 - 0-2)' 



723 



C 



0. 



723 



(3.11) 



This peculiar situation is also encountered in cubic crystals with strong anisotropy, 
such as nickel [25]. In short, the subsonic two-partial surface wave must be slower 
than any in-plane bulk wave (such as the one propagating with speed a/cs/p), but 
is indifferent to the anti-plane wave propagating with speed \fc^p. This principal 
two-partial surface wave is singular because it exists only in the direction Q = 90°, 
although "pseudo-surface waves" may be found in its neighborhood. 

Once X = pv"^ is known, the attenuation coefficients qi are computed from (3.3) 



and the depth profiles follow naturally from (2.16). Figure 3 displays the imaginary 
part of the g^, indicative of the penetration depth, as a function of 9. The horizontal 
straight top (dashed) line is for gi = i. The two other curves (dotted and solid) are 
for the imaginary parts of g2 and q?^. At ^ = 0, there are only two partial modes. 



11 



9(g) 




10 20 30 40 50 60 70 80 90 

Figure 3: Imaginary part of the attenuation coefficients for a surface wave in a 
deformed Mooney-Rivlin material. 

one corresponding to gi = i, the other corresponding to qs = 0.119i. At ^ > 0, a 
third partial mode appears, corresponding to a q2 of the form g2 = i/32, (^2 ^ 0.52. 
In the approximate ranges (0° — 23°) and (68° — 90°), the attenuation factors are of 
the form gi = i, g2 = ^(^2, Is = i/^s with (32 > 0, (3^ > 0. In the approximate range 
(23° — 68°), they are of the form qi = i, q2 = a + i(3, q^ = —a + i/3 with (3 > 0. 
As 9 approaches 90°, and X approaches C5, the imaginary part of the attenuation 
factor gs tends to zero, indicating a deeply penetrating quasi-bulk surface wave 
|26j . Finally at ^ = 90°, a singular two-partial principal surface wave exists, with 
one mode corresponding to gi = i and the other corresponding to g2 = 0.212i (this 
latter value corresponds to a discontinuity in the representation of g2 as a function 
of 9 and cannot be represented on the graph) . 

4 General case 



4.1 Explicit secular equation for surface waves 

In the general case, where no factorization of the propagation condition occurs, a 
different treatment is required. The "method of the polarization vector" , intro- 
duced by Currie |27], refined by Taziev [20], and recently revisited by Ting [28] . 
proves to be a most effective mean for deriving the secular equation as a polynomial 



12 



in X = pv'^ . It relies on the equations 



C/(0) ■ K(")t/(0) = 0, 



(4.1) 



where K^'^' is the symmetric 3x3 lower left submatrix of the A/"" matrix and n is 
an integer. Hence K^^^ = N^ + XI, K^^) = K^^^Ni + NJK^^\ etc. 

Also, for a wave propagating in the symmetry plane of a monoclinic crystal. 
Ting [28] showed that U{0) is of the form 



C/(0) = f/i(0)[l,ia2,/3i 



(4.2) 



where 0^2, Pi are real numbers. We checked that C/(0) is also of this form in 
the present case of a surface wave propagating in a principal plane of a deformed 
material. 



Computing N ^ and N^, we find that K 



(n) 
12 

It follows that the equations (4.1) written for n 
homogeneous system 



K 



(n) 



0, for n 



-1,1,3. 



-1,1,3 reduce to the non- 



-"13 
-"13 

■(3) 



K 



13 



-"33 

-"33 

.(3) 



K. 



33 



-"22 




"2/3i" 




-"11 


-"22 




/3f 


= 


-"11 


-"22 




«i 
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By Cramer's rule, we find 2/3i = Ai/A, /3^ = A2/A where 



(4.3) 



A 



(-1) 

13 


-"33 


-"22 














-(1) 
^13 


-"33 


-"22 
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-(3) 
'■13 


-"33 


-"22 


















-"11 


-"33 


-"22 












Ai = 


-"11 

-"11 


-"33 
-"33 


-"22 
-"22 
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-"13 


-"11 


-"22 










A2 = 


-"13 


-"11 


-"^22 














-"13 


-"11 


-"^22 



(4.4) 



so that 



A? - 4AA2 



0, 



(4.5) 



which is the explicit secular equation for non-principal surface waves in deformed 
incompressible materials. 



13 



K. 



(1) 

22 



X 



-"■11 



Upon inspection of (2.7) and (2.10), we find that K{^ 



'(1) 



K. 



(1) 

33 



X 



/^, 



disposable factor, i^jg , 



X — 7]. Computing N , we find that up to a common 



'(-1) 7^{-l) 



'(-1) 



7^33 , K22 , -f^ii are proportional to polynomials of 



degree 1, 2, 3, 2 in X, respectively. Similarly, computing N , we find that K 



'(3) 
^33 ; 



'(3) 
'■22 ) 



'(3) 
■11 



(3) 
13 ' 



are polynomials of degree 1, 2, 1, 2, respectively. We conclude from 



the definitions (4.4) of A, Ai, A2 that the secular equation (4.5) is a polynomial of 
degree 12 in X = pv^, just as for monoclinic crystals in linear anisotropic elasticity 
|19j . It is too long to reproduce here but it was obtained in a formal manner with 
Maple and with Mathematica. 



The numerical resolution of the polynomial (4.5) yields a priori 12 roots for X 



From these, we discard at once the complex roots, the negative real roots, and the 
roots corresponding to supersonic surface waves (faster than bulk waves). Out of 
the remaining roots, at most one will yield attenuation coefficients gi, q2, qs (with 



a positive imaginary part) from the propagation condition (2.13), such that the 
exact secular equation ( |2.21[ ) is satisfied. 



To conclude this Section, we check that the rationalized secular equation (4.5) 



is consistent with the known secular equation for principal surface waves. When 
6* = 0, it is easy to see that K]^^ 



(-1) 



K 



(1) 



13 



K 



(3) 



13 



0, therefore A = A2 = 0. 
Thus, (4.5) reduces to Ai = 0. By substituting X = 712 — 72iC^) we find that Ai 
factorizes into the product of a quadratic in ("^ and two cubics in (, one of which 



is indeed Dowaikh and Ogden's [5j equation (3.10) 



4.2 Example: Varga materials 

The standard Varga strain energy function [291 EH] is defined as 



Vr = C(Ai + A2 + A3-3), 



(4.6) 



where the material parameter C is constant. This strain energy function has been 
introduced to describe natural rubber vulcanizates. It leads to the following ex- 
pressions for the quantities jij and f3ij, defined in (2.8), 



7i, = CA^/(Ai + A,), A, = C\iXj/i\i + \j). 



(4.7) 



Simple shear plays an important role in the experimental determination of a 
strain energy function [22]. Consider now a half-space made of Varga material, 
subject to an amount of shear 7 along X3, see Figure 4. With a view to the 
possible non-destructive evaluation of sheared rubber, we are interested in the 
propagation of a surface wave in any direction in the plane of shear OXiX^. We 
assume that this plane is free of normal load, thus o"2 = 0. 
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Figure 4: An elastic material under simple shear; the dashed lines represent the 
body at rest. 



The principal stretches associated with the amount of shear 7 may be deter- 
mined using the following equations fl\, 



Ai-Al 



-1 



7: 



Ai>l, 



A, 



A3 = K'x 



2 



(4. 



We substitute the parameters (4.7) into the secular equation (4.5), and then solve 



that equation numerically for the amounts of shear: 7 = 0.5, 1.0, 1.5, corresponding 
to the angles of shear: tan~^7 = 26.56°, 45°, 56.31°, respectively. We obtain the 
velocity of the surface wave propagating in a direction making an angle 6 + ijj with 
CXi, where ip = {l/2)tan-\2/-f) = 37.98°, 31.72°, 26.56°, respectively. Note that 
here OXi is not a principal axis, for the two principal axes Oxi and Oxs in the 
plane of shear are at an angle ip with OXi and with OX^, respectively. 

For the selection of the relevant root of the secular equation out of the 12 
possible ones, we followed the checking procedure described at the end of the 
previous Section. For example, at 7 = 1 and 6 + ip = 0, the rationalized secular 
equation (4.5) has only 2 real roots, giving ^ pv^ jC = 0.848 or ^/pv^JC = 0.884. 
Both roots yield three attenuation factors with a positive imaginary part from 



the propagation condition (2.13). However, the exact secular equation (2.21) is 



satisfied only with the first root. Hence, with a 32 digit precision. Maple finds that 



\nuji/[ujui{m - ujii)] - 1| < 10"^^ when ^/pxpJC = 0.848, indicating that (|2.21[) is 



satisfied; on the other hand, |raci;i/[a;ni(^ — i^ii)] — 1| > 1.96 when ^J pv'^ jC = 0.884, 
indicating that (2.21) is not satisfied then. 

Figure 5 displays the dependence of \J pv^ jC on the angle 6* + V^ over the range 
[0° — 180°] . The solid, dot, and dash-dot curves correspond to an amount of shear of 
0.5, 1.0, 1.5, respectively. The figure confirms what is to be expected intuitively: as 
the half-space is more and more sheared, the strain-induced anisotropy increases. 
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Figure 5: Scaled surface wave speed as a function oi 9 + ip for a Varga material 
subjected to simple shear deformations. 

and its influence on the surface wave speed is more and more marked. It also shows 
that the surface wave travels at its fastest (slowest) speed along the direction of 
greatest (least) stretch, thus allowing for an acoustic determination of the directions 
of the principal stretches in sheared rubber. 
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